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Abstract 

Yang's pure space equations (C.N. Yang, Phys. Rev. Lett. 33, 445 (1974)) gen- 
eralize Einstein's gravitational equations, while coming from gauge theory. We study 
these equations from a number of vantage points: summarizing the work done previ- 
ously, comparing them with the Einstein equations and investigating their properties. In 
particular, the initial value problem is discussed and a number of results are presented 
for these equations with common energy-momentum tensors. 
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1 Introduction 



Over 20 years ago C. N. Yang introduced a system of equations which, while gener- 
alizing Einstein's vacuum equations, take the Yang-Mills equations of gauge theory as 
their underlying structure. These equations, which we refer to as Yang's equations, as 
obvious candidates (at the classical level at least) for the unification of gravitational and 
gauge theory, have been studied from a number of perspectives during the subsequent 
years. Our intention in this paper is threefold: to summarize the work that has been 
done on these equations as they relate to classical general relativity, to argue that these 
equations are worthy of closer scrutiny, and, finally, to extend the work that has already 
been done. 

The basic variable in Yang's theory is, as in general relativity, a Lorentz metric g ab 
on a 4- manifold M. The equations that this metric is required to satisfy are 

V a R bc - V b R ac = 0, (1) 

where R ab is the Ricci tensor and V is the covariant derivative associated with g ab . 

Yang's equations arise naturally by applying the Yang-Mills condition *D*F = to 
the curvature 2-form (Riemann tensor) of an so(3, 1) valued Levi-Civita connection. In 
tensorial language, this reads V a R a bcd = 0, which by the Bianchi identities, is equivalent 
to i). 

We do not propose to allow these equations to supercede those of Einstein as the field 
equations of space-time, but rather intend on the one hand to investigate the consequences 
of Yang's equations as being of interest in their own right, and on the other to study the 
implications of Yang's equation for classical general relativity. 

As a third order non-linear system, Yang's equations are significantly more complex 
than those of either general relativity or Yang-Mills theory, both of which are second 
order. Nonetheless, shortly after Yang published his paper a number of exact solutions 
were discovered by Pavelle @ [§, Thompson |][|, Ni §] and Chen-Long et al §g. 
Later work by Aragone and Restuccia investigated the linearized version of Yang's 
equations, while Mielke (l(J has looked at the equations with emphasis on the double 
anti-self dual ansatz, which yields a gravitational analog of the instanton pseudoparticle 
of Yang-Mills theory. 
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More recently, these equations have been subsumed into a general quadratic la- 
grangian approach, where the basic variables are the metric and the connection (which 
is not necessarily torsion free - see Baekler, Heyl and Mielke [O], Maluf [[12]] Jl3[] and 



Szczyrba [14|). Further interest in the equations has been generated by suggestions that 
they may be suitable for describing the propagation of gravitational waves in a non- 
vacuum (see van Putten |15[|). 

In the next section we shall outline a number of general properties of solutions to 
Yang's equations. Many of these properties have been described in one form or another 
over the last 20 years. In addition, we study the relationship between the curvature 
orthogonality condition (a consequence of Yang's equations) and the algebraic structure 
of the Ricci tensor. 

Any classical dynamical equation is required to be deterministic, in the sense that, 
given some initial (possibly constrained) data on a spacelike slice of spacetime, the equa- 
tions have a solution, at least for a short time into the future. Section 3 discusses this 
initial value problem for Yang's equations and indicates how this can be shown to be 
well-posed. 

Section 4 looks at the compatability of different energy-momentum tensors with 
Yang's equations. In particular, we show that all Einstein-Maxwell fields obeying Yang's 
equations lie in Kundt's class and all perfect fluids have Robertson- Walker geometry. 

There are many open questions concerning these equations which the authors con- 
sider worth pursuing. Describe the asymptotic fall-off for YPS. Does a peeling theorem 
apply? Is there a hamiltonian formulation for the dynamics of the theory? Are the 
equations linearisation stable? How do gauge symmetries of the Levi-Civita connection 
relate to Killing vectors of the metric? Further study of the existing exact solutions is 
also warranted. 

Throughout, we will refer to Lorentz 4-manifolds (M,g a b) which satisfy Yang's equa- 
tion as Yang Pure Spaces (YPS)Q. 



3 The abbreviation shall be applied to both the singular and plural cases. 



3 



2 Some Basic Properties of Yang Pure Spaces 

In this section, we review three properties of solutions of Yang's equation. These prop- 
erties have been noted previously by other authors (see e.g. Thompson @]||); our aim 
is to present (and extend) these results as a preliminary step prior to obtaining solutions 
of Yang's equations under various assumptions. 

Firstly, we note that a contraction of Yang's equations @) yields 

V a R% - d b R = 0, 

so that when this is compared with the contracted Bianchi identity 

Va-^b — \dbR = 0, 

we find that a necessary condition for (M, g a b) be a YPS is that it be of constant scalar 
curvature. 

On the other hand, because in four dimensions the full (second) Bianchi identities are 
equivalent to 

V a C a bcd = V[ c R d ] b - 
Yang's equations (|l|) is equivalent to 

y a C a bcd = 0, R = constant. (2) 

From this, we can immediately identify three large classes of solutions: 

Theorem 2.1 The following are sufficient (though not necessary) for (M,g a b) to be a 
YPS: 

(a) (M,g a b) is an Einstein space, 

(b) (M,g a b) is conformally flat with constant Ricci scalar, 

(c) [M,g a b) has parallel Ricci tensor. 

It is worth noting that Yang's equations are not conformally invariant. Indeed 

Theorem 2.2 Suppose that (M,g a b) is a YPS. Then a conformally related metric also 
solves Yang's equations if and only if the conformal curvature of (M,g a b) is degenerate 
of type N, with twistfree repeated principal null direction. 
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It is also worth mentioning that all solutions of class (b) are obtained by solving the 
nonlinear wave equation 

n M n+ ±m 3 = o, (3) 

for the conformal factor 0, of the metric tensor g a f, = Q 2 rj a b, where O m is the d'Alembertian 
of flat space-time and R is the constant Ricci scalar of g a i>. In fact, the existence of a 
solution to this equation constitutes a version of the "Lorentz Yamabe problem" : 

Lorentz Yamabe Conjecture Let (M,g a b) be a closed or asymptotically flat Lorentz 
manifold, then equation (^) has a solution. 

Note also that because of the Ricci and Weyl tensor terms in the Bianchi iden- 
tities uncouple. Hence in the Newman-Penrose (NP) formalism, Yang's equations are 
equivalent to the contracted Bianchi identities (equations 7.69 - 7.71 of Kramer et al [17] 



or equations 4.12.40-41 of Penrose and Rindler |1{|]), and the remaining Bianchi identi- 
ties (equations 7.61-68 of []17]] or equations 4.12.36-39 of |18[|) with either all Weyl tensor 
terms ignored and A = R/24 held constant, or all Ricci tensor terms ignored. One or 
other version may be used according as one has information regarding the Weyl or Ricci 
tensor. 

Secondly, Yang's equations imply the equation 

V DD '*ABCD = 0, 

for the Weyl spinor. Hence the classic Goldberg-Sachs theorem holds: 

Theorem 2.3 In a YPS (M, g a &) a null vector k a is a repeated principle null direction 
(pnd) of the Weyl tensor if and onlt if k a is tangent to a shearfree null geodesic. 

Thus this theorem, which is an indispensible tool in the study of exact solutions of 
Einstein's equation, may also be used in the study of YPS. 

Thirdly, on taking the covariant derivative of Yang's equations and applying the Ricci 
identities, one can obtain the equation 

which has the irreducible form 

C a b[cdP b e] = 0, ( 4 ) 
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where P a t = R a t — (R/4)g a i, is the trace-free Ricci tensor. Thompson [1] has pointed 
out some of the consequences of this equation, which we will refer to as the Curvature 
Orthogonality Condition (COC), for the Petrov type of non-Einstein YPS. We wish to 
focus on these consequences in more detail. 

In the section 4, we shall investigate the compatibility of different types of physical 
energy-momentum tensor with Yang's equation. It is therefore of interest to determine 
what may be deduced from (Q) , assuming a particular structure for the energy-momentum 
tensor, which determines, via Einstein's equation, the structure of the Ricci tensor. In 
particular, we will determine the Petrov types allowed by the COC for the different Segre 
types of the Ricci tensor. In order to do so, it will be convenient to write (||) in NP form. 
There are nine independent equations as we can see from the following. 

Considering the tensor 

T a bcd — C ae [b c P ^ , 

there are four choices for each of the index sets a and bed. From the symmetries of C a b c d 
and P a b, we find that T a i, c d is completely trace- free yielding the six equations T a acc i = 0, 
and in addition obeys 

± [abc]d 3 A dabci 

which is one further equation. There 8X6 no further possible identities, so that T^cd 

has 

nine independent components in all. Transvecting with appropriate members of a null 
tetrad, these nine equations may be written 



- $21 *0 


+ 2$ 2 2*1 + $02*1 - $0l(*2 + $2) + $00*3 


= 0, 


(5) 


$02*0 " 


- $20*0 + 2$i *i - 2$ i*i - $oo(*2 - *2) 


= 0, 


(6) 


$21*1 


- $12*1 + 2$ii(^ 2 - *2) + $01*3 - $10*3 


= 0, 


(7) 


$22*1 


- $12(2^2 + *2) + $11*3 + $02*3 - $10*4 


= 0, 


(8) 


$22 (*2 


- ^2) + 2*21*3 - 2*12*3 + $02*4 - $20*4 


= 0, 


(9) 


$22*0 " 


- 2$i2*l + $Q2(*2 - *2) + 2$01*3 - $00*4 


= 0. 


(10) 



Notice that equations (|5|), (||), fli~0| ) are complex, giving two real equations each. 

The standard algebraic classification of the Ricci tensor (Hall |ll| , Kramer et al ]17| , 
§5.1) relies upon the number and multiplicities of the eigenvalues, and the nature (space- 
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like, time-like or null) and multiplicities of the eigenvectors of the eigen-problem 

R a bV b = \ v a . 

In order to use (||) in the form (|5|) to (|i0[), we will construct a null tetrad which is based 
on these eigenvectors and consider the NP components of the Ricci and Weyl tensors 
on this tetrad. A sample calculation is given below, and Tables I- IV give a complete 
description of the Petrov types allowed by the COC for a given Segre type. In the 
tables, the first column indicates the Segre type of the Ricci tensor. The second column 
gives the eigenvectors of the Ricci tensor, with time-like, null or complex eigenvectors 
before the comma, and eigenvectors with the same eigenvalue grouped in parentheses. 
Relations which arise for the NP Weyl scalars are given in the third column, and from 
these relations, the allowed Petrov types are determined and indicated in the fourth 
column. Where the Segre type has an important physical interpretation, this is given in 
the final column. Note that the COC provides no information about conformally flat or 
Einstein spaces (Segre type {(1,111)}), and so these will be omitted from the discussion 
below. 

A Ricci tensor of Segre type Al non-degenerate ({1,111} in Segre notation) has a 
pseudo-orthonormal tetrad of eigenvectors, {u a , u„} a= 2,3,4 with corresponding eigenval- 
ues pi,p a , no two of which may be equal. (The pseudo-orthonormality conditions are 
9abVa v a = —9abU a u b = 1, no sum over a.) From these we construct the null tetrad 



k a 



u a + vf) (u a -vf) a _ (v% + ivf) _ a _ (*ug - ivf t 



n„. = = — , m 



s/2 ' a V2 ' V2 ' y/2 ' 

Then we can write 

Rab = ~(pl + P2)k( a n b ) + \(p2 - pi)(k a k b + n a n b ) 

+(p 3 + P4)m( a ' ffl b) + |(P3 - P4)(m a m b + m a fh b ). 

In this case, the NP components $00 = ^22 an d $02 £ R- must be non-zero, and 
may also be non-zero; all other terms are necessarily zero. Specializing equations (|5|) - 
( Jl0| ) to this case, we find immediately that ^1 = ^3 = 0. Furthermore, in order to avoid 
degenerecy among the eigenvalues, we must have {^o = ^2} C R. Then the only 
allowed Petrov types are I and D (c/. §4.4 of Kramer et al [0]). The rest of the table for 
Class Al Ricci tensors is deduced by allowing degenerecies among the eigenvalues above. 
The same form for the Ricci tensor obtains. The results are summarised in Table I. 
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For a Ricci tensor of class A2, there exists a complex conjugate pair of eigenvectors 
z± = k a ±in a with eigenvalues pi±ip2,P2 / 0, and a pair of space-like eigenvectors v%, V4 
as above. k a ,n a are necessarily null, and the normalization k a n a = — 1 may be imposed. 
Then the Ricci tensor may be written as 

Rab = -1p\\a n b) + P2{k a h - n a n b ) + (p 3 + j o 4 )m (a m b ) + ±(p 3 - p 4 )(m a m 6 + rh a fh b ), 

and the following Table II results. No physical Ricci tensor may have this structure. 

For a Ricci tensor of class A3, there exists a double null eigenvector k a with eigenvalue 
pi and a pair of spacelike eigenvectors v$ , v 4 as above. Taking n a to complete a null tetrad 
in the standard way, the Ricci tensor may be written 

Rab = -2pify a n 6 ) + \k a k b + (p 3 + p 4 )m {a fh b) + ±(p 3 - p 4 )(m a m 6 + fh a fh b ), 

where we must have A / 0, for otherwise n a would be a fourth independent eigenvector. 
The results are in Table III. 

A class B Ricci tensor has a triple null eigenvector k a with eigenvalue p\ and a unique 
spacelike eigenvector v\ orthogonal to k a . We complete a null tetrad with n a and v%. 
Then the Ricci tensor may be written 

R ab = -2pik^ a n b) + 2ak( a v 3a) + (p 1 + p 4 )m (a m b ) + ±(pi - p 4 )(m a m 6 + m a m b ), 

where u must be non-zero. Then Table IV results. There are no physical Ricci tensors 
of this type. 

Several results may be read off Tables I-IV. For instance we have the following theo- 
rem. 

Theorem 2.4 If k a is a double (respectively triple) null eigenvector of R ab , then k a is 
at least a 2-fold (respectively 3-fold) repeated pnd of the Weyl tensor. 

Note that the converse of this result is not implied by the COC. For example, for class 
A2, the Weyl tensor is type D if and only if (following the results above) ^0 = ^4 = 0, 
so that k a and n a are both 2-fold repeated pnd's. However the Ricci tensor has no null 
eigenvectors. 

The COC is implied by Yang's equation, but is obviously true for more general space- 
times. In fact we shall see in section 4 that the condition V[ a R b ] c = imposes quite 
severe constraints on space-times having certain Segre types. 
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3 The Initial Value Problem 



The Cauchy initial value problem for the Einstein equations of general relativity has 
received much attention over the last 40 years, from the original work of Lichnerowitz 
p0|| and Choquet-Bruhat PHI , to the sharper results of Hughes, Kato and Marsden p2fl 
and the recent work of Choquet-Bruhat and York p3|]. Similarly, the initial value problem 



for the Yang-Mills equations on Minkowski space has been treated by Kerner [24], Eardley 



and Moncrief p5|][pq] and Klainerman and Machedon [27]. Both the Einstein and Yang- 
Mills equations initial value problems have been solved in somewhat restricted situations. 
For the Einstein equations, short time existence only has been established, except in 
certain cases. Indeed, in view of the singularity theorems (Hawking and Ellis ]^]), it 
appears that long term existence is not possible in the general case. On the other hand, 
while long term existence has been established for the Yang-Mills equations, this has only 
been possible when the gauge group is compact and the base manifold is conformally 
related to Minkowski space-time (Choquet-Bruhat, Paneitz and Segal (^9|). In the case 
of non-compact gauge group, blow-up can occur (Yang and it is not known what 
conditions guarantee long term existence over a general non-conformally flat Lorentz 
manifold. 

We do not expect Yang's equations to have a long term well-posed initial value prob- 
lem for a number of reasons. In the first case, they generalize the Einstein vacuum 
equations, and so the singularity theorems of general relativity will also apply to them. 
Secondly, when viewed from a gauge-theoretic perspective, the equations lie naturally in 
the special orthonormal frame bundle of the base manifold, and therefore the gauge group 
is SO(3,l), which is non-compact. Consequently, we consider only short time existence 
for these equations. 

We assume that the spacetime is the product of V 3 x 3ft, where V s is a compact 
oriented 3-manifold and that the curves {p} x 3ft are timelike everywhere with respect 
to g, while Vt = V s x {t} are spacelike everywhere for all t. Thus if we denote the unit 
normal to Vt by rj then 5 a6 r/ a % = — 1 and g restricted to Vt is positive definite. Endow 
M with a smooth positive definite background metric g and denote the associated Levi- 
Civita connection by V. Let h be g restricted to Vq. Define the Sobolev spaces E S (V 3 x /) 
to be the space of tensors <fi on V 3 xl such that 
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(i) the restriction of cj) and its derivatives V of any order \z\ <s to each Vt is almost 
everywhere defined and square integrable in the metric h. Set 




(ii) the map I— > Sft with ti— ► ||0||.H a _ 1 (Vt) is continuous and bounded, while I— » 5ft with 
ti — >• ||0||j3" s (y t ) i s measurable and essentially bounded. 

E^(V 3 x I) endowed with the norm 

U\\e 3 (V^xI) = Ess suptei \\4>\\ Hs (v t ), 

is a Banach space. For s > §, £' s (y 3 x/) has nice embedding and multiplication properties 
(see Choquet-Bruhat et al pl|| ). Moreover, these are the appropriate spaces for looking 
at the initial value problem for quasi-linear wave equations, such as the Einstein vacuum 
and Yang- Mills equations. We shall use these spaces to study the initial value problem 
for Yang's equations. 

We now find the constraints of Yang's equations on Vo . Denote by 5 g the pullback of 
g to Vq. The second fundamental form K is a symmetric 2-tensor on Vo given by 

Kab = Va Vb^cVd- 

Since Yang's equations are third order in g^, we also require certain 2nd derivatives of 
g a b to be prescribed on the initial hypersurface. Introduce the symmetric 2-tensor Nab 
on Vq which is the pullback of the Ricci tensor of g a b to Vq\ 

Ned = Yc VdRab- 



Proposition 3.1 The equations ([Jj when pulled back to Vq read: 

D [c N b]a = K a[b D d K c]d -K a[b D c] K (11) 
Kf a N b]d = D [a D d K b]d (12) 
constant = 3 R + K 2 - K ab K ab - 2N% . (13) 

Here K = g ab K ab is the mean curvature of the embedding and D is the Levi-Civita 
connection and 3 i? the scalar curvature associated with 3 g. 
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Now we state a theorem which establishes short time existence for Yang's equations. 



The details can be found in Guilfoyle [32|. 

Theorem 3.2 (Local Existence for the Yang's Equations) Let V 3 be a compact 
orientable 3-manifold and 3 g G H^(V 3 ), K G H^iV 3 ), N G H2CV 3 ) be symmetric 2- 
tensors on V 3 , with 3 g positive definite. Let 3 g, K and N satisfy the constraints ({TJj, 
( fJl| ) and Then there exists L C and a unique Lorentz metric g G E^iV 3 x /) 

satisfying Yang's equations with 

g\v = 9, 9 g\v = K and Ricci(g)\ Vo = N. 

Proof: 

The proof follows that of the initial value problem for the Yang-Mills and Einstein 
equations, where the quasi-linear equations in an appropriate gauge are viewed as a 
perturbation of the corresponding linear systems. Then, by use of energy estimates for 
the linear system, a mapping on the spaces E S (V x /) is shown to be a contraction for 
small enough interval / (and large enough s). Since these are Banach spaces, such a map 
has a unique fixed point, the solution of the non-linear system. 

For the Yang system of equations, a number of extra complications arise. Yang's 
equations constitute a third order quasi-linear system of partial differential equations for 
g a b and so fall outside of the usual (2nd order) wave equation framework used for the 
initial value problem of general relativity and gauge theory. In order to overcome this 
difficulty we can separate the connection from the metric and propagate them separately 
by wave equations, to which we can apply standard techniques. The key to the success 
of this method is the fact that if the connection and metric are initially compatible, 
they remain so throughout the evolution. This follows from a geometric identity which 
dictates the propogation of torsion. The appropriate gauge condition for these equations 
is a combination of the harmonic co-ordinate condition 

^(ft-rt) = o, 

and a generalized Lorentz gauge for the connection 1-forms A 

V C A $ = 0. 

c(a) 

With this approach, the proof of the initial value problem goes through as in the Einstein 
case. □ 
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4 The Energy- Momentum Tensor of a YPS 

As outlined above, we do not consider replacing Einstein's equation with Yang's, but 
rather consider the latter to be a complementary condition which is to be imposed on 
space-time. We consider YPS to be worthy of study in their own right, but also feel it 
is important to investigate the interaction of the two theories. Taking this point of view, 
we are led naturally to the following question. What are the consequences of Yang's 
equations for the space-times most commonly studied in general relativity? 

In this section we focus on these consequences for space-times satisfying Einstein's 
field equation, 

Rab - \RQab = 8irT ab , (14) 

with an energy-momentum tensor having one of the following three physically significant 
forms: 

(a) Perfect fluid. 

Tab = (M + p)u a u b + pg ab , u a u a = -1. (15) 

(b) Electromagnetic field. 

T ab = F ac F b c - \F cd F cd , F ab = ~F ba . (16) 

(c) Self-interacting Scalar Field. 

8vrT ab = Va4>V b 4> - ±(V c 0V c </> - V(<f>)) gab . (17) 

Combining Yang's equations with Einstein's leads to 

V [a T 6]c = 0. (18) 

We consider the consequences of this equation for each of the three cases above. The 
energy-momentum tensor always obeys the conservation equation V a T a b = 0, and so ( |i~8| ) 
implies T = g ab T a b = constant. 

4.1 Perfect Fluid 

The condition T = constant leads to the equation of state 

p=±/j + c, (19) 
12 



for some constant c, and the conservation equations for any perfect fluid read 



fi + e(n+p) = o 

V a p + pU a + (p + p)il a = 0, 

where a dot indicates covariant differentiation along the fluid flow lines. Then a straight- 
forward calculation gives 

= u a V [c T b]a = (p + p)(3u [b u c] - V [c u b] ). 

Transvecting with u b then yields 

{p + p)ii a = 0, 

and so assuming /i+p / 0, which we shall do henceforth, we obtain u a = 0, so that the 
fluid flow lines are geodesic. The conservation equations then give 

V a p = 0(p + p)u a , 

so that the spatial gradients of both p and p vanish. Furthermore, this allows us to 
calculate 

= V [c T fc]Q = {p + p){\0g a[b u c] + u [b V c] u a + V [c u b] u a ). 
Transvecting with u b then yields 

V c u a = \9h ac , (20) 

where as usual h ab = g ab + u a u b is the metric tensor on the 3-spaces orthogonal to u a . 
Therefore the fluid flow lines are shear-free and twist-free as well as geodesic. It is well 
known (Krasihski |33|| , Ellis that these are necessary and sufficient conditions for a 
perfect fluid filled space-time to have Robertson- Walker geometry. If p + p = 0, then p 
and p are both constant, giving R ab = constant x g ab , an Einstein space. Since ( p0| ) along 
with (19) is also a sufficient condition on (|l||) to ensure equation (|l8|) , we can summarise 



as follows. 

Theorem 4.1 A perfect fluid space-time (M,g) with p + p ^ is a YPS if and only if 
(M,g) is a Robertson- Walker space-time with p = ^p + c for some constant c. 
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This theorem may be rephrased in terms of Segre types. The usual interpretation is that 
a perfect fluid space-time is one for which the Ricci tensor has Segre type {1, (111)}. 
The energy density p and pressure p are determined from the time-like eigenvalue pi 
and space-like eigenvalue p2 by pi = —Att(p + 3p), p2 = 4:ir(p — p). The non-degenerecy 
condition pi ^ p2 is equivalent to p + p 7^ 0. In addition, p and p are usually required to 
satisfy cert ain positivity requirements, the weak or dominant energy conditions (Hawking 
and Ellis p§, §4.3). Clearly the result ab ove does not depend on such conditions, and 



so may be restated as follows. 

Theorem 4.2 A Segre type {1,(111)} space-time (M,g) is a YPS if and only if(M,g) 
is a Robertson- Walker space-time with constant Ricci scalar. 

Comparing with the first table of §3 above, we see that in this case, the COC is 
quite weak. In fact, since Robertson- Walker space-times are conformally flat, there are 
no space-times of Petrov type /, II, D or N in this class. The central relation which 
constrains solutions to be of Robertson- Walker form is (|20|). From this equation, the 
Ricci identities yield conformal flatness and 3-spaces of constant curvature. The same 
occurs if there is a unit spacelike eigenvector n a of the Ricci tensor obeying 

V a n 6 = l(V c n c )(g ab - n a n b ). (21) 

In this case, the 3-spaces orthogonal to n a are time-like, i.e. have Lorentzian signa- 
ture. Following the steps above, it is easliy seen that (^) holds for a Ricci tensor with 
Segre type {(1,11)1} and R = constant, with n a the spacelike eigenvector with distinct 



eigenvalue. Thus we have the following counterpart of Theorem 4.2 



Theorem 4.3 A space-time (M,g) with Segre type {(1,11)1} is a YPS if and only if 
(M, g) has constant Ricci scalar and line element 

with k constant. 

Finally, we note that the presence of a cosmological constant in Einstein's equation does 
not alter any of the results derived above. The only difference is that the constant c in 
(1£) above is in the general case equal to (4A — R)/3, where A is the cosmological constant 



and R is the Ricci scalar. 
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4.2 Electromagnetic Fields 

We turn next to the case of a YPS generated by an electromagnetic field. The two 
cases, where the electromagnetic field is non-null and where it is null will be dealt with 
separately. We use the NP formalism throughout this section. 

4.2.1 Non-null Electromagnetic Fields 

A non-null Einstein-Maxwell space-time is known as a Rainich geometry and there are 
two conditions which are necessary and sufficient for space-time to be such. These are 
the algebraic part, 

RabRK = \9acRbdR hd + 0, (22) 

and the analytic part, 

V [a a b] =0, a b = (RcdR^r^bgefR'hVfR^. (23) 

The first of these is equivalent to the Segre type of the Ricci tensor being {(1, 1) (11)}, 
while the second follows immediately from Yang's equation. Equation (|2^) allows for 
the existence of a cosmological constant in fll4|), so that while the trace of the energy- 
momentum tensor necessarily vanishes, the Ricci scalar need not. However since we are 
using the original form (14) of Einstein's equation, we wish to rule out this possibility, 



and so demand for the moment that R = 0. We have then 

Theorem 4.4 All Segre type {(1, 1) (11)} YPS with R = are non-null Einstein- Maxwell 
space-times. 

We can explicitly determine all such space-times. In this case, there exists a null tetrad 
field for which 

Tab = 4^i^i{/c (a n fc ) + m (a m b) }. (24) 

k a and n a are the pnd's of the electromagnetic field. Thus the only non-vanishing Ricci 
tensor term is 

$11 = 87T0101 > 0. 



The field equations are the separated Bianchi identities. From these, we immediately 
obtain 

K = a = \ = iy = p = fi = TT = T = 0, 
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^0 = ^1 = ^2 = ^3 = ^4 = 0, 

-D$n = A3>n = <5$n = 0, 

so that space-time is conformally flat (and hence by (Q) includes all conformally flat 
Einstein-Maxwell space-times), $n is constant and both pnd's of the electromagnetic 
field are non-diverging. Hence the space-time lies in Kundt's class (Kramer et al [|l7j 
§27). There is only one such space-time, the conformally flat Bertotti- Robinson solution 
(Kramer et al [17] §10.3), for which the line element may be written 

ds 2 = -Idudr - 2<S> u r 2 du 2 + 2(1 + $nCC)" 2 ^CrfC (25) 

The coordinates x a = {u, r, C, are respectively a label for the surfaces orthogonal to 
k a (k a = V a u), an affine parameter along the integral curves of k a and holomorhpic 



coordinates on the 2-spaces of constant curvature given by u, r constant. In fact (24) 
is the direct product of two 2-spaces of the same constant curvature, and provides an 
example of a decomposable space-time as discussed by Thompson || . This space-time is 
static; £ a = (2<j?nr 2 , 1,0,0) is a hypersurface-orthogonal timelike Killing vector field. 
We can also completely determine all YPS generated, via Einstein's equation with a 



cosmological constant, by a non-null electromagnetic field. Here, (14) is replaced by 

Rab - \RQab + ^9ab = 8TlT ab , 

where A is the cosmological constant. Then 

R = 4A = constant, 



and with T ab as in ( 24 ) , the only difference for the spin coefficients and curvature tensor 
is 

R A 

^2 = = • 

12 3 

Again, the field equations may be completely integrated, yielding a type D (in general) 
space-time which is decomposable into two 2-spaces of constant curvature. The line 
element may be written as 

ds 2 = -2dudr - 2($n - ^)r 2 du 2 + 2(1 + ($ n + ^)(()- 2 d(d(. (26) 



The difference between this and the conformally flat (25) lies in the different values of 
the constant cuvature of the two 2-spaces. 

We can summarise the results for non-null electromagnetic fields as follows. 
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Theorem 4.5 All Segre type {(1, 1) (11)} YPS have line element given by $2q ) and have 
R = 4A. The space-time is generated, via Einstein's field equations with cosmological 
term A, by a non-null electromagnetic field. 

As with any Rainich geometry, the electromagnetic field is only determined up to 



a constant duality rotation, <fii — ► e ta (f)i, a constant. (26) may be generated by an 



electrostatic, magnetostatic or general static electromagnetic field. 

4.2.2 Null Electromagnetic Fields and Pure Radiation 

By definition, an energy-momentum tensor which is one of these types has the form 

T a b = Hk a kf>, 

for some null vector field k a and function PL. Then on any null tetrad based on k a , we 
have from Einstien's equation $22 = 8nH, and all other Ricci tensor terms vanish. (In 
the presence of a cosmological term, R = constant may also be non-zero.) Equivalently, 
the Ricci tensor has Segre type {(2, 11)} with vanishing Ricci scalar. 

From the results of the section 2, we see that the Weyl tensor must be algebraically 
special with repeated pnd k a , which by the Goldberg-Sachs theorem must be shear- free 
and geodesic. The only other consequence of the COC is that ^2 = ^2 for any null 
tetrad preserving the direction of k a . With k = a = 0, Lorentz transformations of the 
null tetrad may be used to set e = p, — p = and r = n = a + j3. Then the only field 
equations (Bianchi identities with Weyl tensor terms ignored) which are not identically 
satisfied are 

5<&22 = -T$22, £>$22 = = p$ 2 2, 

and so p = 0. Thus k a is non-diverging, and the space-time lies in Kundt's class (Kramer 
et al (TtJ §27). For the case of a null electromagnetic field, 4>2 (t ne only non-zero electro- 
magnetic field tensor term) must satisfy the relevant Maxwell equations. These results 
also follow when the cosmological constant is non-zero. Thus we have the following. 

Theorem 4.6 A YPS generated by a pure radiation energy-momentum tensor or a null 
electromagnetic field, or equivalently one with a Ricci tensor of Segre type {(2,11)}, lies 
in Kundt's class. The non-diverging geodesic repeated pnd is the double null eigenvector 
of the Ricci tensor, and in a null tetrad adapted to this vector, ^ = ^2 • 
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The converse of this result ("Kundt's class with null electromagnetic filed or pure 
radiation energy-momentum tensor is a YPS") is not true, as Yang's equation gives 
an 'extra' Bianchi identity arising from the splitting of the identities into homogeneous 
equation for the Ricci and Weyl tensors. Modulo the implications of this extra equation, 



the results of §27.6 of Kramer et al 17] apply to this class of YPS. It may be possible 
to use this extra equation to obtain further results on thes e solutions. This is currently 
being investigated. 

4.3 Self-interacting Scalar Fields 

The energy momentum tensor for a self-interacting scalar field (p with interaction potential 
V(4>) is given by 

87rT ab = <f> a <t> b - g ab (l<Pc(f> c - V (</>)), (27) 

where (ft a = V a (p and for convenience, the gravitational constant has been absorbed into 
the definition of 4> and V. The governing equation for (ft is equivalent to the vanishing 
divergence of (p7|); 

U<j> + V'{4>) = 0- (28) 
Then Einstein's field equation (14) gives 



Rab = 4>a4>b ~ Vg ab , (29) 

and hence 

(j) c c = R + AV. (30) 

Yang's equations imply that R is constant. 

If </> c <^ c = over some open subset of space-time, then we see from ( p0| ) that V = —R/4 
is constant, and consequently the Ricci tensor has Segre type {(2, 11)}, as <fi a is null. This 
situation was dealt with in the previous subsection, and so we assume henceforth that 
(f) c (f) c is non-zero and we can therefore restrict our attention to open subsets of space-time 
on which the sign of ^ a (/> a does not change. We thus note immediately that Theorem 



Q| and Theorem |43| apply in this case; for (j) c (f) c < 0, the Segre type is {1, (111)}, while 
for 4> c (j) c > 0, the Segre type is {(1,11)1}. We now proceed to see how these geometries 
arise. 
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Following the steps of §4.1, we find that the covariant derivative of 4> a is given by 

Va^ = 3(0 b /r V<^ c - V'g ac . (31) 

This equation is equivalent to (|l^) for an energy-momentum tensor of the form ( p7| ) . The 
unit vector field normal to the surfaces E^ of constant is 

n a = (±<p c <f> c r 1/2 r, 

where here and throughout, the upper sign corresponds to 4> c (j) c > (E^ time-like), and 
the lower sign to (p c (p c < (E^ space-like). We find that ([H]) leads to 

V a n 6 = -V'^^T^hab, (32) 

where = g a ^ =F ^a^b is the metric tensor of E^. Using ([H]) and (|32]), an explicit 
expression may be obtained for the second covariant derivative of n a . Then the Ricci 
identities may be used to prove that space-time is conformally flat and that E^ have 
constant curvature. In addition, we find that V must obey 

(R + W)V" - 3(V'f = l(R + 3V){R + W). (33) 

Equation ( p2D allows us to introduce 'co-moving coordinates' (the case </> c </> c may be 
referred to as a tachyon fluid) such that 

n a = ±V a u, u = x°. 

Surfaces of constant u are surfaces of constant <j>, and so 4> = <t>{ u ), V = V{u). The 
possible line elements are 

ds 2 = ±du 2 + A 2 {u) (l + j(x 2 + y 2 =f z 2 ))' 2 (dx 2 + dy 2 T dz 2 ) , (34) 

and k may be normalized to -1, or +1. The 3-spaces E^ are Lorentzian for 4> c <ft c > 0, 
and Riemannian for <fi c (j) c < 0. 

Imposing the condition R = constant on ( |34| ) yields a differential equation for A(u) 
which may be integrated to give 

A 2 (u) = -6/ciT 1 + a exp (J±§ «^ + oi exp f-y / ±f , (R ^ 0) (35) 
A 2 (u) = ±ku 2 + a u + a 1 , (R = 0). (36) 
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In ( |3q ) a coordinate transformation of the form (u — > u+ constant has been used to 
set a constant of integration to zero. (|34| ) along with ( |35|) with the lower sign gives the 
possbible line elements for a YPS with a perfect fluid energy-momentum tensor. 

The scalar field 4> an d potential V may be determined as follows. Using a dot to 
indicate differentiation with respect to u, we have 



b a = ±4>n a , 



and so using (53), 



In the coordinates of (133), we find 



and so (|28|) gives 



□0 = ±( < ^ + 3^~ 1 i(/)), 



which has the first integral 



where a is a constant. Thus by 



+ 2A~ l A4> = 0, 



<p = aA 2 (u), 



V = ±\A- 



(37) 



Then evaluating the terms in (|33|), we find that A must satisfy 



A = -\a 2 A~ 3 ±§A. 



(38) 



We compare this equation with (j35|) to determine a. 

For i? = 0, we find that a = ao- (|37|) may be integrated and we find 



In 



u 



±ku + ao 

For R 0, comparison of (|3^) and (j35l) yields 



V = ±J(±k U 2 + a u)- 2 . 



a 2 = ±12R~ 1 T %a Q axR. 



is determined via (37) and (^) by an elliptic integral. 
Summarising the main results, we have the following. 
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Theorem 4.7 A YPS generated by a self-interacting scalar field (j) obeying V ' a (f)SJ a (j) 7^ 
has line element 

ds 2 = edu 2 + A 2 (u)(l + \{x 2 + y 2 - ez 2 ))- 2 (dx 2 + dy 2 - edz 2 ), 

where e = sgn(V a <pV a 4>) and A 2 {u) is given by $BJ. 



5 Conclusions 

Perhaps the most important result reported here is that the short-time initial value 
probelm for Yang's equations is well-posed. This is a prerequisite for YPS to be worthy 



of study as physical models, and has not been demonstrated before [32]. Other authors 
have emphasised their criticisms of Yang's equations as the fundamental gravitational 
field equations Q, |]35| , This has been done principally on the basis that Birkhoff's 
theorem does not apply to YPS; source-free (V[ a i4]c = 0) spherically symmetric solutions 
are not necessarily static, and indeed the static spherical source-free solutions form a 
four parameter family, rather than the one parameter (Schwarzschild mass) family of 
general relativity. Consequently, one can easily produce 'vacuum' spherically symmetric 
solutions, which are candidates for the gravitational field of the sun, but which display 
unphysical characteristics - lack of gravitational red-shift, lack of bending of starlight, 
and incorrect values (and direction!) for the perihelion shift. 

We have attempted here to establish an alternative philosophy for Yang's equations: 
they are to be used in conjunction with Einstein's equation and/or appropriate boundary 
conditions in our description of space-time. In this way, well established results of rela- 
tivistic astrophysics and cosmology may be maintained. Our hope is that a study of these 
equations may shed further light on the classical gauge theoretic structure of gravity and 
in turn classical gauge theory in general. For example the link between a symmetry of 
a YPS as a gauge configuration and as a space-time may yield insight into the question 
of how the former is to be defined p6| ]. This work is to be undertaken in the future; we 
are currently investigating the structure of static YPS. Details on the short-time initial 
value problem will appear elsewhere. 
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Table I: Class Al Ricci Tensor 



Segre Notation 


Eigenvectors 


Weyl Scalars 


Allowed 
Petrov Types 


Physical 
Interpretation 


{1,111} 


{U, V 2 V 3 Vi} 


*1 = ^3 = 0, 
{tf = *4,*2} CR 


I or D 




{(1,1)11} 


{{u,V 2 )v 3 V±} 


*1 = ^3 = 0, 
{* ,*4,*2} CR 


I II D or N 




{1,1(11)1 


{u,V 2 (v 3 l> 4 )} 


*1 = * 3 = 0, 


I or D 




{(1,1)(H)} 


{(U, V 2 )(V 3 V4)} 


*1 = * 3 = 0, 

^ 2 gR 


I II D or N 


Non-null 
em field. 


{1,(1H)} 


{it, (v 2 v 3 v 4 )} 


#2 G R 


I II D or N 


Perfect fluid; 
scalar field. 


{ (1,11)1} 


{ (u, v 2 v 3 ) V A } 


G R, A =0 - 4. 


I II D III or N 


Scalar field. 



Table II: CLASS A2 Ricci Tensor 



Segre Notation 


Eigenvectors 


Weyl Scalar s 


Allowed Petrov Types 


{zz,U } 


{Z + Z^,V 3 V4} 


= ^3 = 

{^ 4 = -^ ,^ 2 } cR 


I or D 


{zz,(ll)} 


{ z+z-, (v 3 v 4 )} 


^ = ^3 = 
{# 4 = -* ,* 2 } GR 


I or D 



Table III CLASS A3 Ricci Tensor 



Segre Notation 


Eigenvectors 


Weyl Scalar s 


Allowed 
Petrov Types 


Physical 
Interpretation. 


{2,11} 


{k,v 3 v 4 } 


*o = ^1 = *3 = 0, 
{*2,^4} CR 


II D or N 




{2,(11)} 


{k, (v 3 v 4 )} 


*o = ^1 = ^3 = 0, 

# 2 gR 


II D or N 




{(2,1)1} 


{{k,v 3 )v 4 } 


*o = = 0, 

{^2,^3,^4} CR 


II D III or N 




{(2,H)} 


{(k,v 3 v 4 )} 


^ = % = 
^2 G R 


II D III or N 


Null em field, 
pure radiation. 



Table IV: CLASS B Ricci Tensor 



Segre Notation 


Eigenvectors 


Weyl Scalar s 


Allowed Petrov Types 


{3,1} 


{k,v 4 } 


*o = *1 = ^2 = 0, 
{* 3 ,*4} CR 


III or N 


{(3,1)} 


{(k,v 4 )} 


* = *1 = ^2 = 0, 
^3 GR 


III or N 



